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ABSTRACT 

The origin of the different spectral components present in the high-energy (UV to X-rays/gamma-rays) spectra of Seyfert 
galaxies is still being debated a lot. One of the major limitations, in this respect, is the lack of really simultaneous 
^ ' broad-band observations that allow us to disentangle the behavior of each component and to better constrain their 

interconnections. The simultaneous UV to X-rays/gamma rays data obtained during the multiwavelength campaign on 
(f) . the bright Seyfert 1 Mrk 509 are used in this paper and tested against physically motivated broad band models. 

■ Mrk 509 was observed by XMM-Newton and INTEGRAL in October/November 2009, with one observation every four 

' days for a total of ten observations. Each observation has been fitted with a realistic thermal Comptonization model 

for the continuum emission. Prompted by the correlation between the UV and soft X-ray flux, we used a thermal 
Comptonization component for the soft X-ray excess. We also included a warm absorber and a reflection component, as 
required by the precise studies previously done by our consortium. The UV to X-ray /gamma-ray emission of Mrk 509 
can be well fitted by these components. The presence of a relatively hard high-energy spectrum points to the existence 
of a hot (kT ~ 100 keV), optically-thin (r ~ 0.5) corona producing the primary continuum. In contrast, the soft X-ray 
component requires a warm (kT ~ 1 keV), optically-thick (r ~ 10-20) plasma. 

Estimates of the amplification ratio for this warm plasma support a configuration relatively close to the "theoretical" 
configuration of a slab corona above a passive disk. An interesting consequence is the weak luminosity-dependence of 
its emission, which is a possible explanation of the roughly constant spectral shape of the soft X-ray excess seen in 
AGNs. The temperature (~ 3 eV) and flux of the soft-photon field entering and cooling the warm plasma suggests that 
it covers the accretion disk down to a transition radius Rt r of 10-20 R g . This plasma could be the warm upper layer of 
the accretion disk. 

In contrast, the hot corona has a more photon-starved geometry. The high temperature (~ 100 eV) of the soft-photon 
field entering and cooling it favors a localization of the hot corona in the inner flow. This soft-photon field could be 
part of the comptonized emission produced by the warm plasma. In this framework, the change in the geometry (i.e. 
Rtr) could explain most of the observed flux and spectral variability. 

Key words, galaxies: individual: Mrk 509 - galaxies: active - galaxies: Seyfert - X-rays: galaxies 



1. Introduction 

Since the discovery of a high-energy cut-off in the X-ray 
emission of Seyfert galaxies, first detected by SIGMA and 



OSSE in NGC 4151 (|.Iourdain et all . Il992t iMaisack et all 
I1993D an d then commonly observed in about 50% of Seyfert 
galax ie s (IZdziarski et al.l Il993t iGondek et al.l. 19961 : Matt, 
120011 : iPerola et all 120021: iDeluit fc Courvoisierl [2003; 
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iDadinal . 120071: iMolina et aLl . l2009ft . the primary continuum 
of this class of AGN is commonly believed to be thermal 
Comptonization of soft seed photons produced by a 'cold 
phase" , presumably the accretion disk, on hot (~ 100 keV) 
thermal electrons (the "hot phase" also called corona). 
The different Compton scattering orders produced through 
this process eventually add t ogether to produce a roughl y 
cut-off power-law shape (e.g. ISunvaev fe Titarchuk| [r980'). 
which agrees well with the high-energy (>2 keV) spectra 
of Seyfert galaxies. 

The nature and origin of this thermal hot plasma are, 
however, still largely unknown, mainly owing to the lack 
of sensitive instruments at high energies. This hinders our 
ability to put tight constraints on the spectral shape and 
to disentangle the primary continuum from the reflection 
component. At lower energies (< 2 keV), the presence of 
deep absorption features, the signatures of ionized plasma 
(the so-called warm absorber, WA hereafter) surrounding 
the central engine, and soft X-ray emission in excess of 
the extrapolation of the high-energy continuum to these 
energies (the so-called soft X-ray excess) does not allow a 
direct view of the primary emission. Thus, we are generally 
left with the 2-10 keV band, which is far too narrow for 
a deep understanding of the origin of the high-energy 
broad-band emission. 

The spectral and flux variability, close to the time scale 
of a day for black-hole mass of 10 8 M Q but covering a wide 
frequency range, makes things even more confused if we are 
not able to separate the contribution of the different spec- 
tral components. For example, the X-ray spectral shape of 
the continuum is expected to harden (the photon index T 
decreases) when the corona te mperature (hence th e high- 
energy cut-off) increases (e.g. lHaardt et al.l I1997D . Such 
variations in the X-ray spectral shape may be produced 
by intrinsic changes in the corona properties (for example, 
changes in the heating process efficiency) and/or by varia- 
tion in the environment like changes of the soft-photon flux 
(and consequently of the corona cooling) produced by, e.g., 
the accretion disk. 

The origin of the soft X-ray excess has been strongly 
debated. It was rapidly proposed that it could result from 
warm Compton izatio n in an optically thick plasma (e.g. 
IWalter fe Finklll993HMagdziarz et aLlll998l . and see more 
recently iDone et al. M2012D . In this case, a correlation be- 
tween the UV and soft X-ray emission is expected. If it 
is p roduced by blurred io nized reflection, as proposed by, 
e.g.. lCrummv et al~l (|2006ft . then the soft X-rays should be 
linked with the hard (> 10 keV) X-rays and reflection com- 
ponents; i.e., a broad iron line and a bump at ~ 30 keV, are 
expected. Other authors have ex plained the soft X-ray ex - 
cess by smeared absorption fe.g. iGierlinski fe Donei [2006'). 
but the observations apparen tly require blurring eff ects that 
are too large to be realistic ([Schurch et all [2009) . Finally, 
the effect of partial covering ionized absorption can also 
account for part of the observed soft excess without re- 
quiring material to be moving at extreme velocities (e.g . 
iMiller et al1l2008l and the review of iTurner fc Milled l2009j) . 

In summary, it is obvious that improving our under- 
standing of Seyfert galaxy high-energy emission requires 
long and broad-band monitoring. 



Mrk 509 is a Seyfert 1 galaxy that harb ors a super mas- 
sive black hole of 1.4 x 10 8 solar masses (jPeterson et all 
2004). It was the object of an intense optical/UV/X- 
rays/gamma-rays campaign in late 2009, with the par- 
ticipation of five satellites (XMM-Newton, INTEGRAL, 
Swift, Chandra, and HST) and ground-based telescopes 
(WHT and PAIRITEL). The sour ce characterist i cs an d 
the whole campaign are explained in lKaastra et al.l (|2011[ ). 
The main objective of this campaign was to study in detail 
the physical properties and structure of the warm-absorber 
outflows in Mrk 509. The heart of the monitoring are ten 
simultaneous XMM-Newton/TNTEGRAL observations, 
with one observation every four days (60 ks with XMM- 
Newton and 120 ks with INTEGRAL). This paper is part 
of a series focusing on different results obtained with the 
campaign. For example, the study of the stacked 600 ks 
RGS spectrum enabled us to put strong constraints on 
the outflow nature and geometry dDetmers et all 120111 : 
iKriss et all . 120111: iKaastra et all . 12012ft . 

Cl osely related to the present paper, Mehdipour et al. 
(|201lL hereafter Mil) have studied the evolution of the 
UV/X-ray broad-band emission of Mrk 509 using the OM 
and pn instruments. Mil show that the UV and soft X-ray 
continuum emission are tightly correlated, suggesting 
a common origin for the two continuum emission com- 
ponents. This favors an origin through Comptonization 
instead of blurred absorption or emission, as discussed 
above. This interpretation has recentl y been suggested 
for several narrow-line Seyfert galaxies (iDewanean et all 
[20071: [Middieton et al.l . [200alJin et alll2009ft . The UV /soft 
X-ray correlation observed in Mrk 509 agrees with the 
soft X-ray emission being the upper-energy part of the 
Comptonization of the UV photons in a warm (kT ~ 1 
keV) and optically thick (r ~ 10-20) corona. It is worth 
noting that these results do not rule out the presence 
of other components in the soft X-rays (such as blurred 
reflection or absorptions). But they do not dominate the 
variability on a time scale of the few days observed during 
this monitoring. 

The study done by Mil was limited to data below 
10 keV, and the primary continuum was supposed to 
be a simple power law. It is known, however, that real 
Comptonization spectral shapes may be si g nificantly 
differ e nt from a pow er law (e.g. lHaardt! [1991: [S vensson 
119961 : iPetrucci et al.l [2000 i. Moreover, the interpretation 
of the spectral variability in terms of physical quantities, 
such as the temperature and optical depth of the corona, 
is not straightforward from the photon index and high 
energy cut-off of a phenomenological cut-off power-law 
approximation, so it requires the use of approximate 
formulae to estimate these parameters. 

In the present paper we aim at testing realistic 
Compt onizat ion models in a similar way to the Petrucci 
et al. (|2004l) analysis of the simultaneous IUE/RXTE 
campaign of NGC 7469. The main advantages of the Mrk 
509 campaign are that the data cover the X-ray range 
from the soft X-ray bands up to 200 keV, while the RXTE 
energy range, in the case of NGC 7469, was limited to 
3-20 keV. Furthermore, the XMM-Newton data allow us to 
put strong constraints on the optical/UV continuum and 
the iron line complex and to analyze the WA structure 
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Fig. 1. Count rates light curves in different energy ranges. 
From top to bottom: 2120 A, 0.5-1 keV, 3-10 keV, 20-60 
keV, and 60-200 keV. The horizontal solid lines are the 
average values and the horizontal dashed lines represent 
the corresponding ± la uncertainties. 



with Swift, started on September 4, 2009, for a total of 19 
pointings, and the last observation, with Chandra, ended 
on December 13, 2009. The HST/COS data were taken 
about 20 days after the end of the simultaneous XMM- 
Newton/IN TEGRAL observatio ns and consist of ten HST 
orbits (see iKaastra et al.l 1201 ll for a detailed description 
of the campaign). The archival FUSE data were obtained 
October 15, 2000, i.e. about nine years before our campaign. 



We refer the reader to iPonti et al.l (|2012ft (P12 
hereafter) for details about the XMM-Newton/ pn data re- 
duction and especially the gain calibration and correction. 
This data reduction was done with the latest version of the 
SAS software (version 10.0.0), starting from the ODF files. 
The contribution from soft proton flares was negligible 
during the 2009 monitoring. The XMM-Newton/OM data 
treatment and the host galaxy subtraction are detailed 
in Mil, and we also refer the reader to this paper for 
more information. INTEGRAL/ISGRI data were reduced 
with OSA 8, with additional ghost-cleaning improvements 
provided in OSA 9. The presence of SPI annealing meant 
that 40% of the observations cannot be used for the SPI 
instrument. This unfortunately results in too low a signal- 
to-noise ratio for the SPI data to be used. Observation 6 
was split in two owing to operational contingency, however 
we added both parts to improve statistics. 

The total XMM-Newton/ pn spectrum and the pro- 
duction of the corresponding response and ancillary files 
were performed with the MATHPHA, ADDRMF, and 
ADDARF tools within the HEASOFT package (version 
6.8). The total INTEGRAL spectrum, as well as the 
corresponding response and ancillary files, were generated 
using the OSA spejpick tool. The total on-source exposure 
time is ~ 430 ks for XMM-Newton/pn and ~ 620 ks for 
INTEGRAL. 



precisely. Thus, for the first time, we are able to do a 
broad-band spectral analysis of UV to X-ray data of a 
Seyfert galaxy, on a month's time scale, including all the 
major spectral components known to be present in this 
class of objects. 

The paper is organized as follows. We briefly describe 
the observation and data reduction in Sect. [21 referring 
the reader to the other papers of the collaboration when 
appropriate. In Sect. [31 we present light curves in differ- 
ent energy bands and interband correlations including the 
INTEGRAL data. In Sect. [31 a principal component anal- 
ysis of the UV/X-ray/7-ray data is shown. The detailed 
spectral analysis is done in Sect. [5j The results are dis- 
cussed in Sect. [Hand we state the conclusions in Sect. [7] 

2. Observations and data reduction 

For details of the observations of XMM-Newton, 
INTE GRAL, and Swift, precisely please see IKaastra et all 
(pOllT) . The XMM-Newton/INTEGRXL monitoring 
started on October 18, 2009 and ended on November 20, 
2009. During this campaign, we also obtained UVOT, as 
well as Swift X-ray, Chandra, and HST/COS data. We 
further more use archival FUSE data. The first observation, 



3. Light curves and flux correlation 

The XMM-Newton and INTEGRAL count rate light curves 
are reported in Fig. [T]for different energy ranges: at 2120 
A i.e., the wavelength of the OM UVW2 filter, 0.5-1 keV 
and 3-10 keV for XMM-Newton/ pn, and 20-60 keV and 60- 
200 keV using INTEGRAL ISGRI data. A peak-to-peak 
variability of 30-40% is visible in the XMM-Newton/pn light 
curves and is also detected in the UV band (see also Mil). 
Some variability seems to be present in the 20-60 keV and 
60-200 keV band of ISGRI. Noticeably, the INTEGRAL 
20-60 keV light curve looks similar to the XMM-Newton 
3-10 keV one. However, given the error bars, the ISGRI 
soft- and hard-band count rate light curves are formally 
consistent with being constant. 

An important result observed during this campaign is 
the strong correlation between the UV and soft X-ray (<0.5 
keV) flux, while no correlation was found between the UV 
and the hard (>3 keV) X-rays (see Mil). This suggests 
that the UV and soft X-ray excess variability observed on 
the timescale of the campaign is produced by the same 
spectral component, possibly thermal Comptonization (see 
below). 
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Fig. 2. Count rate-count rate correlation coefficients be- 
tween different energy bands for the 10 observations. In 
each panel, the green star is the XMM-Newton/OM light 
curve in the 3-6 eV band; the three blue crosses correspond 
to the XMM-Newton/pn soft X-ray bands: 0.2-0.5 keV, 0.5- 
1 keV, 1-2 keV; the three red diamonds correspond to the 
XMM-Newton/pn medium X-ray bands: 2-4 keV, 4-6 keV, 
6-12 keV; and the two black triangles correspond to the 
INTEGRAL /ISGRI hard X-ray bands: 20-60 keV and 60- 
200 keV. In each panel, the light curves are correlated with 
the one marked with the vertical dotted lines. We have 
overplotted the 95% and 99% confidence levels for the sig- 
nificance of the correlation (it corresponds to correlation co- 
efficients of 0.64 and 0.78, respectively, for 10 points) with 
horizontal dot-dashed and dashed lines, respectively. 



We have plotted in Fig. [5] the linear Pearson correlation 
coefficients between light curves in different energy ranges, 
including those above 10 keV from the INTEGRAL/ISGRI 
instrument. The error bars reported in this figure were es- 
timated in the following way. To simplify, we focus on the 
error of the Pearson correlation coefficient of 0.78 obtained 
between the 3-6 eV and 0.2-0.5 keV count rate light curves. 
For both energy ranges, we simulated 200 light curves of ten 
count rates, each of these count rates being drawn from a 
Gaussian distribution probability centered on the observed 
values (in these energy ranges) and with a standard de- 
viation equal to the 68% count rate error of the observed 
data. 

We take the example of the simulation of the faked 3-6 
eV count rates for Obs 1. The observed count rate for Obs 
1 is 14.06 cts s _1 with a 90% error of ±0.76 cts s -1 . Thus 
we draw 200 faked 3-6 eV count rates from a Gaussian 
centered on 14.06 cts s _1 and with a standard deviation of 
0.76/1.65=0.46 cts s _1 . We proceed in the same manner 
for the 0.2-0.5 keV energy band. Then we correlate the 
200 simulated 3-6 eV light curves with the 200 simulated 
0.2-0.5 keV light curves to obtain a set of 200x200 fake 
correlation coefficients. We define the correlation coefficient 
Cmin (respectively c max ) as the value below (respectively 
above) which we have 5% of the simulated correlation 
coefficients. In the present example, c m i n = 0.47 and 
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Fig. 3. The first (left panels) and second (right panels) 
principal components of variability. The upper panels show 
the spectra corresponding to the maximal and minimal 
coordinate cxj t k over the j=10 observations of the first 
(k=l,left) and second (k=2, right) eigenvectors. The mid- 
dle panels show the ratio of the maximum and minimum 
spectra to the total spectrum of the source. The bottom 
panels show the contribution of each component to the to- 
tal variance as a function of energy. 



0.85, and we use c„ 



and 



as the corre- 



sponding 90% error on the observed correlation coefficient. 
We apply the same procedure for the other correlation 
coefficients shown in Fig. [2] We have also overplotted, 
in this figure, the 95% and 99% confidence levels for the 
significance of the correlation. This correspond to correla- 
tion coefficients of 0.64 and 0.78, respectively, for 10 points. 

Figure [2] shows that the UV and soft X-rays (<1 keV) 
bands are indeed well correlated even if the errors on the 
correlation coefficient are relatively large. But, in contrast, 
the absence of correlation with the higher energy bands 
points to separated components between the low (<1 keV) 
and high (>1 keV) energies. Noticeably, the XMM-Newton 
energy bands above ~2 keV appear to vary in unison, and 
also correlate well with the 20-60 keV INTEGRAL energy 
band, suggesting a common physical origin. 



4. Principal component analysis 

Before going deeper into the spectral analysis, we perform 
a principal component analysis (PCA) to search for vari- 
ability patterns. PCA is a powerful tool for multivariate 
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data analysis which is now widely used in astronomy 



Paltani fc Walter] 119961: 



Vaughan & Fabian 



iMalzac et alj|2006t iMiller et alj|2007t iPris et al.l l2011h . It 
allows transforming a number of (potentially) correlated 
variables into fewer uncorrelated variables called principal 
components. In the present case, the long exposures per 
observation enable the use of a large number of energy 
bins without being limited by the Poisson noise. We choose 
a number of 91 energy bins, from the UV (6 bins) up 
to the X-rays (85 bins), with the INTEGRAL data left 
out due to their low statistics. Consequently, we have 
p=10 spectra obtained at different times t±, ti, t p , and 
binned into n=91 energy bins Ei—i >n . Thus we have a 
p x n matrix F whose coefficients are given by the energy 
fluxes of each spectrum in each energy bin. The PCA 
gives the eigenvalues and eigenvectors Vk=i,n(E) of the 
correlation matrix of F, thus defining a new coordinate 
system that best describes the variance in the data. The 
observed spectra F(E,tj) can be decomposed along the 
eigenvectors with a given set of coordinates ak following 



F(Ei,tj 



fc=i 



Oik(tj)Vk(Ei) for each energy bin Ei = 



X.n- 



The coordinates a.), vary from observation to observation 
and their fluctuations account for the observed variance. 
On the contrary, the eigenvectors are constant and define 
the different variability modes of the set of observations. 

The first PCA components are those that reproduce 
most of the observed variability, with the others dominated 
by the statistical noise in the data. We have plotted the 
two first principal components of our analysis in Fig. [3] 
The first component (left panels of Fig. [3]) consists mainly 
in a variability mode dominated by flux variations in the 
hard energy range (> 1 keV). Most of the sample variance 
(>87%) is produced by this component. The second PCA 
component (right panels of Fig. [3]) is dominated by variabil- 
ity below 1 keV, and represents 10% of the total variance. 
The 4% of variance left is distributed between the higher 
order components not shown here. This result nicely sup- 
ports the idea that two main components are dominating 
the broad-band variability of this campaign, one in the UV- 
soft X-rays and the other in the hard X-rays. 

5. Spectral analysis 

Here we study the spectral energy distribution between the 
optical, UV, X-ray, and gamma-ray bands for Mrk 509. The 
XMM-Newton/^ spectra have been rebinned in order to 
have a minimum of 25 counts per bin and three bins per 
resolution element. The cross calibration normalization be- 
tween XMM-Newton and INTEGRAL is fixed to one. All 
spectral fits were performed using the XSPEC software (ver- 
sion 12.3.0). In Sect. 15.11 the reported errors are at the 
90% confidence level for one interesting parameter, while 
in Sect. 15.21 owing to the potential strong correlations be- 
tween the model parameters, the reported errors are at the 
90% confidence level for two interesting parameters. 

5.1. Above 3 keV 

We first concentrate our analysis on the data above 3 keV, 
i.e. on a domain known to be dominated by the primary 
continuum. This agrees with the PCA analysis (Sect. 0J, 
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Fig. 4. Top: Total XMM-Newton (pn, black) and 
INTEGRAL (ISGRI, red) data above 3 keV. The solid 
line is the best-fit power-law model. Bottom: Ratio 
data/model. The iron line is clearly visible and no high- 
energy cut-off is required. 



which suggests that a unique component dominates the 
variability in this energy range. In this section, we are inter- 
ested in the global spectral behavior of the source. We thus 
simply use a cut-off power-law for the continuum. More re- 
alistic Comptonization models are used in Sect. I5~2l when 
the broad-band UV/X-rays/gamma-rays energy range is 
used. 

To start, we analyze the total spectrum over the entire 
campaign to determine the best model components to be 
used. Afterwards we apply the corresponding model to the 
different observations in Sect. 15.1.21 



5.1.1. The total spectrum 

The best fit and data/model ratio of the total spectrum are 
plotted in Fig. 2) It clearly shows the presence of an iron 
line near 6.4 keV. The presence of a high-energy cut-off is 
not required by the data, and a value below ^200 keV is 
even excluded at the 90% confidence. 

The iron line complex observed during our campaign 
is discussed in detail in PI 2. It can be decomposed into 
a narrow core {a ~ 0.03 keV) and a broader component 
(a ~ 0.22 keV). When looking to archival (XMM-Newton, 
Chandra, and Suzaku) data, the narrow core appears to 
be constant on a several-year timescale, peaking at an 
energy consistent with 6.4 keV. Thi s suggests an origin 
in remote neutra l material (see also iDe Rosa et al.l [2004; 
iPonti et al.l [2009) . On the contrary, the intensity of the 
broader component seems to follow the variations of the 
continuum throughout the campaign (i.e. on a timescale 
of a few days), without a measurable time lag (see Fig. 
4 of P12). Given the sampling of our monitoring, this 
suggests an origin within a few light days from the black 
hole, the width of the line being compatible with a 
Keplerian motion at 40-1000 gravitational radii from the 
1.4xl0 8 Mq supermassive black hole of Mrk 509. Note 
that the relatively large EW of ^50 eV of this component 
suggests that part of it could originate in the inner BLR 
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6.58±0.04 


300±50 


2.9±0.4 


55±20 


2.9±0.3 


468/243 


Obs. 1 


1.71±0.04 


>90 


6.60±0.15 


370^° 


3.0±1.0 


eotf 


(f) 


213/251 


Obs. 2 


1.76±0.02 


>310 


6.82±° ™ 




4.4± ;? 


100+.?° 


(f) 


278/251 


Obs. 3 


1.76±0.04 


>300 


6.55±0.10 


210±5*° 


2.5±1.0 


40±30 


(f) 


287/251 


Obs. 4 


1.74±0.03 


>250 


6.70±0.20 


660+^° 


4-8±?;| 


110+?° 


(f) 


270/251 


Obs. 5 


1.75±0.03 


>100 


6.57±0.15 


360±f 


o Q+1.5 
°- z -1.2 


60±|8 


(f) 


264/251 


Obs. 6 


1.73±0.03 


>120 


6.52±0.05 


<140 


1 Q+0-8 

L -°-o.e 




(f) 


300/251 


Obs. 7 


1.77±0.02 


>350 


6.50±0.10 




O Q+l.l 

°- u -1.0 


sots; 


(f) 


310/251 


Obs. 8 


1.76±0.02 


>200 


6.60±0.10 








(f) 


256/251 


Obs. 9 


1 7fi+ 002 

1 -'°-0.04 


>160 


6.51±0.09 


330+™° 






(f) 


281/251 


Obs. 10 


1.68±0.04 


>80 


6.54±0.08 




9 o+l.O 
z -°-0.7 


80±« 


(f) 


226/251 



Table 1. Best-fit parameters when fitting the data above 3 keV. A cut-off power law is used for the primary continuum, 
a Gaussian for the broad-line component, and a neutral reflection (pexmon) for the narrow iron line and the reflection 
bump. For the individual observations, the PEXMON parameters are fixed (f ) to the best-fit values of the total spectrum. 



(|Costantini. E. et al. in preparation! l2012f ). 

Following P12, we fit the iron line profile with two 
model components. Fo r the narrow core we used the pex- 
mon model of xspec (Na ndra et all . [2007h . It reproduces 
the reflection com p onent following the pexrav model 
(|Zdziarski fe CoppI 119911) of XSPEC, but in addition it 
includes the Ka and K/3 neutral iron line, the Compton 
shoulder, and the nickel Ka line, all these components 
being computed "self -consistently" according to the 
I George fc Fab ian (1991) simulations. We fixed the inclina- 
tion angle to 30 deg. in pexmon. Since it is supposed to 
model remote reflection, the illuminating spectrum should 
be an average of the different states of the central source. 
Thus, we chose to fix the illuminating power-law photon 
index to 1.9 and a fixed high-energy cut-off at 300 keV. 
The exact values of these parameters have no significant 
effects on the results shown in this paper. We let only the 
normalization of pexmon free to vary in the fit of the total 
spectrum. Concerning the broader component we used a 
simple Gaussian whose energy, flux and width are left free 
to vary during the fit. 

The best-fit parameter values for the total spectrum 
are reported in the first row of Table [T] The power-law 
photon index T = 1.74 ± 0.01, the broad line component 
peaks at 6.58±0.04, and has an EW of 55±20 eV and the 
PEXMON normalization is N pexmon = (2.9 ± 0.3) x 10~ 3 . 
The PEXMON normalization is mainly constrained by the 
intensity of the core of the narrow iron-line component. The 
fit is not satisfactory (x 2 /dof=468/243), mainly due to the 
residuals at low energies (close to 3 keV) and possibly due 
to the influence of the soft X-ray excess. 

5.1.2. Individual observations 

We applied the same fitting procedure to the ten different 
observations. The only varying parameters are the photon 
index and normalization of the power law and the high- 
energy cut-off and the characteristics (energy, width, and 
intensity) of the broad- line component. The pexmon nor- 
malization is not expected to vary during the campaign (as 
the signature of remote reflection) and is fixed to the best- 
fit value obtained with the total spectrum. 



The best-fit parameters for each pointing are reported in 
TableQ] For the individual pointings, we always obtained an 
acceptable fit. The best-fit photon index is only marginally 
variable with a range in between 1.68 and 1.77, and the 
presence of a high-energy cut-off is never required. The 
shape of the high-energy spectra of Mrk 509 thus changes 
only slightly during the monitoring. The broad-line com- 
ponent is always statistically needed. The line equivalent 
width is consistent with a constant with an average value 
of 55± 10 eV. The line flux follows the continuum flux 
variability, in agreement with the results obtained by P12. 
However, its energy is only marginally consistent with the 
6.4 keV neutral energy, possibly indicating the presence of 
ionized iron (see P12 for a more detailed discussion). 

5.2. From UV to hard X-rays/soft gamma-rays 
Comptonization models 

By fitting above 3 keV we ignore the spectral complexity 
of the soft X-ray range where soft X-ray excess and 
warm-absorber (WA) features are present. These spectral 
features have to be included in the fitting procedure of the 
whole dataset from UV to hard X-rays/soft gamma-rays. 
They are discussed in some detail in the next paragraphs. 

It is worth noting, at this stage, that the purpose of the 
present paper is not to test different models (especially con- 
cerning the soft X-ray excess or the primary continuum). 
We instead make the choice of fiting a few, physically mo- 
tivated spectral components, which will permit us to dis- 
cuss the corresponding physical interpretation of the data 
in depth. 

5.2.1. The model components 

The warm absorber: the WA has been comprehensively an- 
alyzed by iDetmers et al.l (|201ll ). and we refer the reader 
to this article for more details. We produced a table from 
the best model of the WA derived from this analysis and 
included this table in XSPEC. 

The soft X-ray excess: as already explained in the intro- 
duction, different models exist in the literature to inter- 
pret this component. The most commonly used in the last 
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few years are blurred ionized absorption or reflection and 
warm Comptonization. That no evidence for changes in the 
X-ray WA was found in Mrk 509 despite a soft X-ray in- 
tensity increase of ~ 60 % in the middle of our campaign 
(|Kaastra et all . 1201 lh disfavors ionized absorption as the 
origin of the optical-UV/soft X-ray correlation. 

Concerning ionized reflection, i t is known to b e able 
to produce a soft X-ray exce s s (ICrummv et al. 



iPonti et all . 120091: ICerruti et al.1 1201 H iNoda et ~ 



Moreover, since part of the X-ray emission is reprocessed 
in opt-UV in the disk, we also expect some correlation 
between these energy bands. Two arguments weaken, 
however, the importance of ionized reflection in the 
opt-UV/X-ray variability behavior observed in Mrk 509. 
First we expect a reflection bump peaking near 10 keV 
that should vary with the soft part of the reflection. But 
no correlation is observed between the soft and the hard 
X-ray bands (see Fig. It is true, however, that this 
could be partly explained by the low statistics above 10 
keV. But more importantly, part of the iron line should 
also be produced by this reflection component and would 
be expected to vary with the soft X-rays. However, neither 
the narrow core of the line nor its broader component 
show variations that agree with those observed in the soft 
X-rays (see Sect. I5.1.T]) . 

Following Mil, we chose to model the UV-soft X-ray 
emission by thermal Comptonization. We used the nth- 
COMP model of XSPEC. NTHCOMP covers a wide range of 
parameters and is well adapted to warm (kT~ 1 keVfl 
and optically thick (r ~ 10-20) plasma (i.e. what we call a 
warm corona in the following), such as the one expected to 
reproduce the steep soft X-ray emission. The free parame- 
ters of NTHCOMP are the temperature of the corona kT wc 
0, the soft-photon temperature Tj,f, >wc , and the asymptotic 
power-law index of the spectrum. We prefer to use this 
model instead of COMPTT (which is similar to the COMT 
model in SPEX used by Mil), because the seed photons in 
COMPTT only have a black-body distribution. With NTH- 
COMP we can choose a multicolor-disk distribution. It is 
also fast to use in XSPEC (see also App. [Alfor a comparison 
of these two models) . 

Note however that NTHCOMP does not provide the 
corona optical depth directly. It has to be deduced in an- 
other manner, generally by comparing its spectral shape 
to another model that provides the optical depth. This is, 
however, a possible source of uncertainties and is discussed 
in detail in App. [S] In the following, the optical depth of 
the WARM corona is deduced from a comparison with the 
eqpair model. 



The primary continuum: concerning the high-energy con- 
tinuum, the cut-off power-law model is known to be a poor 
approximation of realistic thermal Comptonization models. 
First, the high-energy cut-off shape of Comptonization 
spectra differs from an exponential cut-off. Besides, 



1 In XSPEC, the lower limit of the temperature parameter of 
NTHCOMP is 1 keV. However, this can be modified and put to 
a lower value. We have checked that the model still works for 
kT~0.5 keV by comparing its spectral shape with the one ob- 
tained with EQPAIR. 

2 The letters "wc" means warm corona and refer to all the 
parameters of the NTHCOMP model. 



Comptonization spectra result from the sum of multiple 
bumps, each one corresponding to the different Compton 
scattering orders experienced by the soft seed photons in 
the corona. For small optical depth, the spectral shape 
is bumpy. Finally, Comptonization spectra naturally 
include the seed photons that cross the corona but are not 
Compton-scattered (the so-called zeroth scattering order). 

We used the C ompto nization model COMPPS 
(jPoutanen fc Svenssonl . I1996T ) of XSPEC to model the 
high-energy primary continuum. COMPPS is well adapted 
to reproducing the emission of a hot and optically thin 
plasma, i.e., what we call the hot corona in the following. 
The fit parameters of COMPPS are the temperature of 
the corona fcTft,J3 , either the optical depth r^ c or the 
kT 

parameter y^c = 4 ^t^ c , where ytc is equivalent to 

m e c z 

the Compton parameter for small (< 1) optical depths, 
the temperature of the soft photons produced by the cold 
phase kTbb.hc (also assuming a multicolor-disk distribution) 
and the geometry of the disk-corona configuration. 

We prefer to use the Compton parameter instead of the 
optical depth in the fitting procedure to minimize the model 
degeneracy. Indeed, the temperature and optical depth are 
generally strongly correlated in the fitting procedure, with 
the same 2-10 keV slope obtained wi th different combina - 
tions of these two parameters (see e.g. lPetrucci et al.ll200l[ ). 
We also assume a slab geometry for the corona-disk system. 

The iron line: following the results obtained in Sect. 
15.11 we use PEXMON to fit the narrow core of the iron 
line and a Gaussian to fit a broad-line component if needed. 

It is important to point out that no link between 
the Comptonized spectrum and the soft UV emission is 
imposed a priori. The model simply adjusts its parameters, 
independently of each other, to fit the data. It is only 
a posteriori that the resulting best-fit values can be 
interpreted in a physically motivated scenario. 



5.2.2. The COS-FUSE data set 

The use of non-simultaneous optical-UV to hard X-ray 
data could produce spurious constraints on Comptonization 
models and should be avoided; however, a few arguments 
make the use of the COS and FUSE data in our spectral 
analysis reasonable. First, the UV continuum shape of Mrk 
509 did not change si gnificantly, even fo r a 60% lower flux as 
observed with FOS ijkriss et all 1201 lh . This suggests that 
the UV spectral shape of Mrk 509 only weakly changes as its 
flux varies. Moreover, some of the Swift data were obtained 
simultaneously with the COS observations in a bandpass 
comparable to the OM bands. The Swift fluxes are at values 
comparable to OM points near the beginning and the end 
of the XMM-Newton campaign, and only 12% be low the 
mean flux of the whole campaign (|Kriss et all 1201 il ). Given 
the lack of evidence for variation in the spectral shape and 
the fact that only a small-scale factor is needed to bring 



3 The letters "he" mean hot corona, and refer to all the pa- 
rameters of the compps model. 
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Fig. 5. Unfolded best- fit model of the total data spectrum 
using a thermal Comptonization component (COMPPS in 
xspec) for the primary continuum. The data are the black 
crosses. The solid black line is the best-fit model including 
all absorption (WA and Galactic). The dashed lines corre- 
spond to the different spectral components, hot (in red) 
and WARM (in blue) corona emission as well as the reflec- 
tion (magenta) produced by pexmon including the effects 
of absorptions, while the dotted line are absorption free. 



region of the accretion flow compared to the WARM corona; 
however, a 100 eV temperature is large for a "standard" 
accretion disk around a ~ 10 8 Mq black hole. This point 
is discussed in Sect. 16.2. II 

The value of 100 keV for the HOT corona temperature 
is compatible with the lower limit on the high-energy 
cut-off E c >200 keV obtained when fitting with a cut-off 
power-law shape for the primary continuum (see Sect. I5.l| ). 
Indeed, as said in the introduction, the cut-off spectral 
shape of a Comptonization spectrum differs from an expo- 
nential cut-off and, roughly speaking, E c is generally about 
a factor 2 or 3 larger than the "real" corona temperature. 
Thus E c > 200 keV agrees with a corona temperature of 
-100 keV. 

The fit of the total spectrum is still statistically not 
acceptable, however, with the presence of strong features 
in the data/model ratio, especially in the soft X-rays. The 
FUSE points are also slightly below the best-fit model (see 
Fig. [5]). A possible reason for these discrepancies, given that 
the fits are reasonably good for each individual observation 
(see next section), could be the source spectral variability. 
The total spectrum is then a mix of different spectral states 
that cannot be easily fit with "steady-state" components. 



the COS data in line with the mean OM fluxes, we decided 
to use of the COS data in our fit procedure. 

The scaling factor to match the FUSE data with 
the OM one, still assuming a constant spectral shape, is 
much larger (~ 2, see Mil for more information on the 
method used to rescale the COS and FUSE data) and 
the use of the FUSE data is more arguable. However, 
it allows for a better constraining of the maximum 
of the UV bump. For these different reasons, we chose 
to include the COS and FUSE data in our spectral analysis. 



5.2.3. The total spectrum 

First, we assume the same soft-photon temperatures 
for the two Comptonization models, i.e., Tbb,hc=Tbb,wc- 
The fit is very bad with x 2 /dof =8850/691, and the 
data/model ratio shows several features in the X-rays. 
Relaxing the constraint of equal soft temperatures between 
the two plasmas enables us to reach a much better fit with 
X 2 /dof =1387/690. The corresponding best-fit parameter 
values are reported in the first row of Table We have 
plotted in Fig.[5]the corresponding unfolded best-fit model. 

The two plasmas clearly have different characteristics. 
The warm corona is optically thick with t wc ~20, a tem- 
perature kT wc on the order of 600 eV and a soft-photon 
temperature Tbb, W c ^3 eV. In contrast, the HOT corona is 
optically thin t>, c — 0.6, with a temperature kT^c of about 
100 keV and a soft-photon temperature Tbb,hc~100 eV. 

In this model, the optical/UV data are entirely fitted 
by the low-energy part of the warm plasma emission (see 
Fig. [5]). The soft-photon bump of the hot corona, instead, 
peaks in the soft X-ray range and contributes to the soft 
X-ray excess. The soft-photon temperature of —100 eV 
suggests that the hot corona is localized in the inner 



5.2.4. Individual observations 

We repeated the same fitting procedure for the ten different 
observations, letting Tbb,hc and Tbb, W c free to vary inde- 
pendently. As in Sect. 15.1.21 the pexmon normalization is 
fixed to the best-fit value obtained with the total spectrum. 
The best-fit parameters for each pointing are reported in 
Table [2] The 90% contour plots of the temperature vs. the 
optical depth are plotted in Fig. [5] for both coronae and for 
the ten observations. The errors reported in Table [2] are 
computed from these contour plots and not from the error 
command of xspec. Indeed, the shape of the \ 2 surface 
is relatively complex for the models used, and confidence 
intervals are found better by determining the locus of 
constant A\ 2 than by the use of the error command. 

The time evolutions of the different best-fit parameters 
are plotted in Fig. [Jj For each parameter we have overplot- 
ted their weighted average value with the corresponding 
±1ct uncertainties. All the light curves are inconsistent 
with a constant at more than 99% confidence except for the 
temperature, optical depth, and soft photon temperature 
of the warm corona. 



5.2.5. Heating and cooling 

We call L to t the total corona luminosity and L s the 
intercepted soft-photon luminosity, i.e., the part of the 
soft-photon luminosity emitted by the cold phase that 
effectively enters and cools the corona. The difference 
Ltot — L s then corresponds to the intrinsic heating power 
provided by the corona to comptonize the seed's soft 
photons. The ratio L to t/L s is called the Compton amplifi- 
cation ratio. 
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Obs Tuic kT wc Twc fcThc Dhc ^~hc Tfjb^iuc ^6(>,/ic ^pexmon \ / dof 







(keV) 




(keV) 






(eV) 


(eV) 


(x 10" 3 ) 




Total 


9 EQC+U.UU5 
Z.Oi)O_ 005 


n w+ u - U2 


18.4+1M 


123 + | 


63+ uul 
u - DJ -o.oi 


o.estn; 


o 4 +o.i 
J -^-o.i 


ioo+i 


2.9 


1387/690 


Obs 1 


2 S6+ 003 
^• JD -o.oi 


n ,10+0.07 

0.48_ 04 


20.3l ;l 




63+ 06 
u - DJ -o.oi 


0.89t™ 


2 6+ 02 
z -°-o.i 


97+1 3 


(f) 


396/415 


Obs 2 


o cc+0.01 

z - JJ -o.oi 


U - M -0.06 


i9.3±?:S 


122±Jf 


63+ 02 


66+ 03 


9 y+0.2 

z -'-o.i 


95+ 2 


(f) 


510/415 


Obs 3 


9 cq+0.01 
Z - JO -0 .01 


U - M -0.06 


1 q Q+ 1 - 2 

15 - J -0.4 


120l? 9 


60+ 02 


65+ ' 07 
u - DO -o.oi 


2 9+ 02 
z - y -o.i 


IO8+2 


(f) 


501/415 


Obs 4 


2 60+ 01 


n ^fi+ 05 

U - OD -0.04 


17 1+ 2 - 1 


164+« 


u - 04 -o.oi 


42+ 02 

u - y:Z -0.03 


q 9+0.1 
°- z -0.5 


102+2 


(f) 


500/415 


Obs 5 


9 K7+0.01 
z - ' -0.01 


U - OU -0.03 




141±» 


n ci +0.01 

u - 01 -o.oi 


4fi+ 02 


q q+0.2 
o -°-0.3 


105±i 


(f) 


467/415 


Obs 6 




u - dy -0.08 


18.5±" 


108+2! 


66+ 03 

u - DD -0.04 


n 70+O.O8 


2 9+ 01 


1161a 


(f) 


518/415 


Obs 7 




U - M -0.08 


19.0± ; 3 8 


H9l 2 3 o 


60+ 02 

u - DU -0.03 


67+ 07 

u -°' -0.02 


2 9+ - 4 

z - y -0.2 


104+2 


(f) 


513/415 


Obs 8 


9 cn+0.02 
Z - JZ -0.01 


u - JZ -0.03 


20.0±g'g 




64+ 005 


1 - UO -0.12 


2 6+ 02 
z -°-o.i 


ii4i 2 3 


(f) 


494/415 


Obs 9 


9 rr+0.02 

z - JJ -o.oi 


4Q+0-04 
u -^ y -0.04 
n cq+0.05 


20.2+?'i 




60+ 02 


67+ 08 

u -°' -0.02 


3 0+ 01 

O.U_ Q 2 


100l 2 


(f) 


450/415 


Obs 10 


9 cq+0.01 




102+^2 


67+ 03 

u -°'-0.03 


n 8/1 +0.04 
U -° 4 -0.05 


2 6+ 03 

z -°-0.2 


108+1 


(f) 


409/415 



Table 2. Best-fit parameters of the total spectrum and the different observations using the Comptonization model 
COMPPS for the primary continuum and nthcomp for the UV-soft X-rays. The optical depth t wc of the warm corona is 
estimated from eqpair fits (see App. [5]) . 



From our best fits obtained for each observation, we can 
estimate the seed's soft-photon luminosity L s that enters 
both coronae. Indeed, the number of photons is conserved in 
the Comptonization process. Consequently, from the num- 
ber of photons measured in our nthcomp and COMPPS 
model components and the temperature of the soft-photon 
field, we can deduce (assuming a multicolor-disk distribu- 
tion) the corresponding soft-photon luminosity that crosses 
and cools each corona (see App. [Bj. In this reasoning, we 
assume that no extra component participates in the ob- 
served opt-UV emission apart the Comptonization spectra. 
Any extra component would decrease the effective seed soft- 
photon luminosity L s , implying that our estimates of L s are 
strict upper limits. 

The total corona luminosity L to t can be deduced from 
the observed luminosity L f, s of each Comptonization model 
component and L s by the relation L to t — 2L f, s — L s e~ T 
(see App. [B]), where r is the corona optical depth. This 
relation assumes that the corona emissions are isotropic. 

We have plotted the values of L s and L to t for both 
coronae, as well as the amplification ratios for the ten 
observations, in Fig. [5] The errors on L a b s and L s (and 
consequently on L tot ) were estimated, in the case of Obs 1, 
from a large number of coronae model parameters (kT wc , 
F wc , kThc, Vhc and the model normalizations), taking into 
account their likelihood given by the contour plots shown 
in Fig. |51 The corresponding flux histograms were fitted 
by a Gaussian from which the 90% error is deduced. These 
errors are on the order of 5% for L t, s and L s . Because this 
procedure is time-consuming, we assume the same order of 
magnitude for the errors for all the observations. 

While the total flux emitted by the HOT corona varies 
during the monitoring, with a significant decrease between 
Obs. 3 and 4, its variation is limited to about 25%. This is 
smaller than the variation in its seed soft-photons luminos- 
ity, which reaches a maximum in Obs. 5, 70% higher than 
the luminosity at the beginning or the end of the monitor- 
ing. Concerning the warm corona, L s and L to t vary sim- 
ilarly, also reaching their maximum in Obs. 5, but ^35% 
higher than the luminosities estimated in Obs. 1 and 10. 

The warm and hot coronae have clearly different am- 
plification ratios (see Fig- IHJs) ■ This is remarkably constant 



and equal to ^2 for the warm plasma as expected given 
its large optical depth (see Eq. IB. 141) and in between 5-15 
for the hot corona, with a significant decrease between 
Obs 3 and 6. 



6. Discussion 

6.1. The disk-corona geometry 

As already said in Sect. 15.21 no link between the comp- 
tonised spectrum and the soft UV emission is imposed a 
priori in our fitting procedure. The model simply adjusts 
its parameters to fit the data. We can, however, compare 
our best-fit values to theoretical expectations and constrain 
our own scenario. 

The spectral shape of Comptonization models depends 
on the temperature and optical depth of the corona. 
Another important parameter is the geometry of the disk- 
corona system, which can be related to the amplification 
ratio. In steady states, the corona temperature adjusts, 
for a given optical depth, to keep this ratio to a constant 
value. For instance, in the case of a plane-parallel corona 
geometry in thermal equilibrium above a passive disk, 
this amplification ratio is expected to be on the order 
of 3 for a thin (r < 1) coron a and 2 in the thick case 
(e.g. lHaardt fc Maraschil |1991| and see App. [B]). Then, 
photon-starved (respectively photon-fed) geometries, i.e. 
geometries where the corona is undercooled (respectively 
overcooled) compared to the plane geometry, have am- 
plification ratios higher (respectively lower) than these 
theoretical values. 

We have overplotted in Fig. [5] the predicted relation- 
ships in steady state between the temperature and the op- 
tical depth assuming a soft-photon temperature of 10 eV 
and different Compton amplification ratios. In agreement 
with the estimated amplification ratios plotted in Fig. [HJs, 
the disk-HOT corona system is clearly in a "photon-starved" 
configuration, with amplification ratios above 3. This sug- 
gests a disk-corona geometry where the solid angle under 
which the corona "sees" the accretion disk is small. This 
is a common result for Seyfert galaxies. Indeed, the slab 
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Fig. 6. Contour plots of the temperature versus optical 
depth for the hot (top) and WARM (bottom) corona for 
the 10 observations. We have overplotted the predicted re- 
lationship t vs. kT in steady state assuming a soft-photon 
temperature of 10 eV and for different Compton amplifica- 
tion ratios: dotted line: L tot /L s =2; dashed line: L tot /L s =3; 
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Fig. 7. Time evolution of the different fit parameters. From 
top to bottom: the hot corona temperature kT^c, Compton 
parameter y^ c , optical depth The, and the WARM corona 
temperature kT WCl photon index T WCl optical depth t wc , 
and the soft photon temperatures Tbb tWC and Ttb.hc- The 
solid lines show the best-fit constant values and the solid 
lines the ±lcr uncertainties. 



geometry cannot reproduce hard X-ray spectra (r2_io <2) 
like those observed in Seyfert galaxies because of the strong 
cooling that is expected fr o m the accretion dis k in such ge- 
ometry (e.g. lHaardtj[l993t lHaardt et all 11994) . A photon- 
starved disk-corona configuration means less cooling com- 
pared to the slab geometry hence harder spectra. 

In contrast, the warm corona agrees very well with a 
configuration of an optically thick plasma covering a passive 
disk. This is the direct consequence of a large optical depth 
(r ffic > 1) and our assumption to fit the Opt-UV and soft X- 
ray data with a unique Comptonization component (see Eq. 
IB. 141 with no disk intrinsic emission) . The slight disagree- 
ment, in Fig. [51 of the predicted relationship in steady state 
Ltot/Lg = 2 and the contours of the temperature vs. the 
optical depth comes from the fact that the theoretical pre- 
diction, correct for optically thin media, is slightly incorrect 
at high optical depth. 



6.2. The physical origin of the warm and hot coronae 

The sum of the seed soft-photons luminosities L s<wc + L St h c 
(dominated by L SyWC ) crossing both coronae is on the order 
of ~ 10 45 erg s _1 for the total spectrum. Assuming that it 
is produced by a multicolor accretion disk around a super- 
massive black hole of 1.4x10 8 Mq, with an inner-disk radius 
of 6R g (R g the gravitational radius associated to the black 
hole) and an inclination angle of 30 degrees, its expected 
inner-disk temperature should be on the order of 15-20 eV. 
This is significantly more than the soft-photon temperature 
~3 eV of the WARM corona but significantly less than the 
soft-photon temperature ~100 eV of the hot corona. 

With the same reasoning, the soft-photon temperature 
of ^3 eV of the WARM corona corresponds to a disk radius 
of ^30 R g . Given the expected slab geometry of the disk- 
warm corona system, it suggests that the warm corona 
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Fig. 8. a) Total flux L to tl^Edd emitted by the hot corona 
and b) intercepted soft photon luminosity L s / LEdd entering 
and cooling the HOT corona, c) and d), like a) and b) 
but for the WARM corona, e) Compton amplification ratio 
Ltot I L s for the WARM (bottom curve) and HOT (top curve) 
coronae, respectively. The dotted line in e) is the expected 
value for a thick corona in radiative equilibrium above a 
passive disk. All the luminosities are in Eddington units 
assuming a black- hole mass of 1.4 x 10 8 A/q . The error bars 
on L b s and L s have been estimated as on the order of 5%. 



is preferably located at this distance on the disk. On the 
other hand, the high temperature of the hot corona and 
the temperature of ^100 eV of its seed soft photons suggest 
a hot plasma in the inner region of the accretion flow. 

6.2.1. warm corona: the disk surface layer? 

That we manage to fit the Opt-UV and soft X-ray data with 
a unique Comptonization component of large optical depth 
implies an amplification ratio of the warm corona equal to 
2 (see Eg. IB. 14]) . i.e., the theoretical value expected in the 
case of a thick corona above a passive disk. This strongly 
suggests that the regions of the accretion disk responsible 
for the optical-UV emission, which acts as the source of 
seed photons for the warm plasma, are essentially heated 
by the warm plasma itself. 

With an optical depth of ~ 20 and a temperature of 
500-600 eV, this (plane) WARM corona could be seen as the 
warm skin layer of the accretion disk. In this case the soft 
emission peaking at 3 eV would come from regions deeper 
in the disk and heated by the warm layer. Such vertical 
stratification of the disk is clearly expected in the case of 



irradiated accretion disks, the upper layer being heated 
to the Compton temperature of the illuminating flux (e.g. 
Nayakshin, Kazanas & Kallman 2000; Nayakshin 2000; 
Ballantyne, Ross & Fabian 2001b, Done & Nayakshin 
2001). 

The optical depth of ~ 20 deduced from our fits ap- 
pears, however, to be well above the expected theoretical 
values, on the order of a few units. This could be explained 
by the very simple assumptions of our corona models, 
which assume a constant temperature and optical depth 
for the whole corona. The inclination of the accretion disk 
could also play a role, since it would increase the effective 
optical depth compared to the theoretical estimates, which 
are computed in the vertical direction of the layer. 

A natural origin of the illuminating radiation pro- 
ducing the warm corona could be the hot corona itself 
(see next section). However, the hot corona luminosity 
is about a factor 2 lower than the total warm corona 
luminosity (sec Fig. [5]); so irradiation is only part of the 
heating mechanism of the upper layers of the accretion 
disk forming the warm corona. Besides, the amplification 
ratio of 2 disfavors intrinsic emission produced by the 
disk deeper layers but instead suggests that the warm 
corona has its own heating process. If this interpretation is 
correct, it means that the accretion heating power is most 
likely to be liberated at this disk surface layer and not in 
the disk interior. Interestingly, recent 3-D radiation-MHD 
calculations, when applied to protostellar disks, show that 
the disk surface layer can indeed be heated by dissipating 
the disk magnetic field, while the disk interior is "mag- 
netically dead" and cold er than a standard accretion disk 
(|Hirose fc Turned. l20ll . The same processes could occur 
in the disk of Mrk 509, giving birth to the warm corona 
needed in our fits. 

It is worth noting, however, that these conclusions 
are model-dependent, and for instance, they strongly 
depend on our assumption to fit the optical-UV and soft 
X-ray emission with a unique Comptonization component. 
The addition of a pure multicolor-disk component in the 
optical-UV band would necessarily decrease the impor- 
tance of the warm corona emission in this energy range 
and consequently would decrease its amplification ratio. 
From a statistical point of view, there is no need to add 
such a multicolor disk component, and thus, with this limit 
in mind, we can conclude that the WARM corona likely 
covers and heats a large part of the accretion disk. 



6.2.2. hot corona: the hot inner flow? 



The total power emitted by the hot corona is on the order 
of 1.3xl0 45 erg s _1 i.e. ~5% of the Eddington luminosity 
of a supermassive black hole of 1.4x10 8 Mq. There are a 
few theoretical models that explain the presence of such 
hot thermal plasma at such high luminosity: for in stance 
the lu minous hot accretion flow (LHAF) proposed bv I Yuan! 
1200 ll which is an extensio n of the ADAF so lutions to high 
accretion rate (see als o lYuan et al.l [200 6). Jet emitting 
disks (JED) solutions (jFerreira et all l2006t ) also provide 
a hot and optically thin plasma. Since no powerful and 
collimated jet is observed in Mrk 509, the JED solution 
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Fig. 10. A sketch of the accretion flow geometry in the inner region of the galaxy Mrk 509. The hot corona, producing 
the hard X-ray emission, has a temperature of ~ 100 keV and optical depth ~ 0.5. It is localized in the inner part of 
the accretion flow (R < Ri n ) and illuminates the accretion disk beyond Ri n , helping to form a warm layer at the disk 
surface. This warm component has a temperature of ~ 1 keV and an optical depth ~ 15 and produces the optical- UV 
up to soft X-ray emission. It extends over a large part of the accretion flow, heating the deeper layers and comptonizing 
their optical-UV emission. In return, part of this emission enters and cools the hot corona. 





0.014 




x> 


0.012 




LU 
_l 

U 
-C 


0.01 




_l 


0.008 






0.006 






0.16 0.18 0.2 0.22 0.24 

Ltot,wc/l-Edd 

Fig. 9. Correlation of the seed soft-photon luminosity of the 
hot corona, L Sy hc, and the total luminosity of the warm 
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is apparently not appropriate. Aborted jets have been 
proposed , however, to explain the properties of Seyfert 
galaxies (IHenri fe Petrucci Il997t iPetrucci fe Henri Il997t 



iGhisellini et all . l200lT A JED could then be present but 
without filling all the conditions for producing a "real" jet, 
such as its insufficiently large radial extension. The HOT 
corona could be present but not the jet. 

The 100 eV temperature for the soft seed photons of the 
hot corona is a factor ~5 above the value expected from 
a standard accretion disk. A few physical arguments could 
help to solve this inconsistency. First, the disk spectra are 
certainly more complex than the simple sum of blackbodies 



due to the possible effects of electron scattering in compar- 
ison to absorption. Special and general relativistic effects 
could also play a significant role on the observed emission. 
These effects lead to modified black-body spectra with an 
"observed" temperature higher than the e ffective temper- 
ature by a factor between 2 and 3 fe.g. ICzernv fe Elvisl 
[19871: lRo^s~eTaDll992T: iMerloni et all 12000ft . The presence 
of the warm corona above a large part of the disk can 
also strongly modify the disk emission , as shown in our 
fits (ICzernv fc Elvisi . [l987t ICzernv etHl l2003t IDone et all 
I2012h . 

An interesting interpretation could be that the seed soft 
photons of the hot corona are those emitted by the warm 
plasma. The warm plasma emission peaks at a few eV and 
extends to a few 100's of eV, so it can be "seen" by the 
hot corona as a soft photon field with an intermediate 
temperature of ~ 100 eV. Such an interpretation cannot 
be modeled easily with XSPEC, but it is nicely supported 
by the strong correlation (with a linear Pearson correla- 
tion coefficient equal to 0.92, which corresponds to a >99% 
confidence level for its significance) between the seed soft 
photon luminosity of the hot corona, L s> hc, and the total 
luminosity of the warm plasma L totyWC (see Fig- [HI) - 

6.2.3. A tentative toy picture 

From the above discussions, a tentative toy picture of 
the overall configuration of the very inner regions of Mrk 
509 could be the following: a hot corona (T e ~ 100 keV, 
r=0.5) fills the inner part of the accretion flow, below 
say Ri n . It illuminates the outer part of the disk and 
contributes to the formation of a warm upper layer (T e ~ 
1 keV, r=20) on the disk beyond Ri n , the warm corona. 
This warm corona also possesses its own heating process, 
potentially through dissipation of magnetic field at the 
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disk surface. It heats the deeper layers of the accretion 
disk, which then radiates like a black body (at Tbb ~3 eV). 
This warm corona upscatters these optical-UV photons 
up to the soft X-ray, a small part of them entering and 
cooling the hot corona. The small solid angle under which 
the hot flow sees the warm plasma agrees with the large 
amplification ratio of ~5-10 estimated from our fits. The 
100 eV temperature obtained by our xspec modeling 
would be explained by the fact that XSPEC would try to 
mimic the WARM plasma emission, which covers an energy 
range between ~ 1 eV and ~ 1 keV, with the multicolor- 
disk spectral energy distribution of the seed soft photons 
available in COMPPS. A tentative sketch is shown in Fig. 1101 

Clearly, in this picture the two coronae should be pow- 
ered by the same accretion flow and should consequently 
be energetically coupled. Assuming that the warm corona 
covers the accretion disk between R ou t and Ri n and then 
that the hot corona fills the space between Ri n and Ri SCO , 
then the accretion power liberated in the WARM corona is 



p 

J ac 



GMM 



i 



Rj 



R„ 



(1) 



while the accretion power liberated in the HOT corona is 



p 

J acc. 



h c 



GMM 

^Risco 



Ri. 



Ri. 



so the ratio between the two is 



P 



acc. he 



Rin 1 



Ri. 



1 



Ri, 
Ro* 



(2) 



(3) 



Neglecting the advection, which seems reasonable given the 
accretion rate of the system, the accretion power liberated 
in each corona is close to their total luminosity. The ratio 
Pacc,hcl 'Pacc,wc is then on the order of L tot ,hc/ L to t,wc and 
Eq. [3] gives a one-to-one relationship between the total 
luminosities L to t,hc and L tot>wc (or, equivalently, L to t,hc 
and L to t,hc/ Ltot.wc) and the inner disk radius Ri n and 
total accretion rate M, as soon as Ri SCO and R ou t are given. 

We have mapped, in Fig. [TTJ the Ri n -M plane into 
the Ltot.hc~Ltot.hc/ Ltot.wc plane assuming R lsco =&R g and 
R O ut=W0Rg- The observed total luminosities agree with 
Eq. [1] and Eq. [2] if Ri n is between 8 and 10 R g and the 
accretion rate M between 20 and 30% of the Eddington 
rate (assuming an accretion efficiency r\ — (2R isco / R g )~ r ). 
Rin strongly depends on Ri SCO and decreases as Ri SC o 
decreases. A value of 10i? g can then be seen as a rough 
upper limit. On the other hand, M only slightly depends 
on Risco- Both Ri n and M very slightly depend on R ou t 
unless R out become very small (< 10-20 R g ). 



6.3. A luminosity-independent spectral shape for the soft 
X-ray excess 

That the warm corona, which is at the origin of the Opt- 
UV to soft X-ray emission in our modeling of Mkr 509, has 
naturally a constant amplification ratio has an interesting 
consequence. It implies a spectral emission that is almost in- 
dependent of the AGN luminosity. Indeed, the photon index 




0.2 0.4 0.6 0.8 

i-tot,hc /Lot ,wc 



Fig. 11. Contour plots of the inner disk radius Ri n (solid 
lines, in units of R g ) and the total accretion rate M (dashed 
lines, in Eddington units and assuming an accretion ef- 
ficiency r\ = (2R lsco /R g Y 1 ) in the Ltotjic Ltot,hc/ Ltot,wc 
plane (see Eqs.[T]and[5]). The points indicate the positions 
of the 10 observations. 



of the comptonized spectrum is directly linked to the ampl i- 
fication ratio (e.g. iBeloboro dov 1999; Malz ac~et al.ll200lh . 
A constant amplification ratio then implies a roughly con- 
stant spectral index. This also depends on the soft photon 
temperature, but this effect is weak as far as this tempera- 
ture is not varying too much. Moreover, in radiative equi- 
librium the corona temperature and optical depth follow 
a univocal relationship. For a given optical depth r corre- 
sponds a unique temperature kT and, to the extent that 
the radiative equilibrium is insured, there is no reason for 
r and kT to vary. The temperature and optical depth of the 
warm corona are indeed consistent with a constant during 
the monitoring (see Fig. [7]). Thus, in these conditions, we 
expect an almost constant spectral shape for the soft X-ray 
excess, regardless of the luminosity of the AGN, in very 



good agreement with observations (e.g . Gierlinski fc Dong 



good agreem ent with observati ons le.g. l Uier 
120041: ICrummv et all 1200^ iPonti et afll200fift 



6.4. Geometry variations 

Simulations by Malzac & Jourdain (|200ll ) show that the 
time scale for a disk-corona system to reach equilibrium, 
after rapid radiative perturbation in one of the two phases, 
is of a few corona light crossing times i.e. shorter than a 
day. We can thus (reasonably) assume that each of the ten 
spectra of Mrk 509 correspond to a disk-corona system 
where hot and cold phases are in radiative balance with 
each other. In this case, a fixed disk-corona configuration 
in radiative equilibrium should correspond to a constant 
heating/cooling ratio. As already said, this is what we 
observe in the case of the WARM corona, its amplification 
ratio being necessarily (given our model) equal to 2 (see 
Fig. |8^). In the case of the hot corona, however, the 
amplification ratio decreases in the middle of the campaign 
(Obs 4 and 5) suggesting a variation in the disk-HOT 
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corona configuration during the monitoring. 

A decrease in the amplification ratio means that the 
disk-HOT corona geometry becomes less photon-starved. 
During Obs 4 and 5 the total luminosity of the WARM 
corona increases (cf. Fig. [5]:), while the total luminosity of 
the hot corona decreases (cf. Fig. |5Ji). If the above picture 
is correct, these variations should be linked somehow to 
variations in the inner accretion disk radius Ri n and/or 
of the total accretion rate. Looking to Fig. [TTJ variations 
of ~5-10% for these parameters are enough to explain the 
observed luminosity variations. The variation in Ri n can 
also explain the variation in the hot corona amplification 
factor. 



6.5. A pair-dominated hot corona? 

From the luminosities L to t and L s , we can estimate 
the compactness ltot = Ltot&T / '(Rm e c 3 ) and l s = 
L s ut I (Rm e c 3 ) where ctt is the Thomson cross section, m e 
the electron mass, c the speed of light, and R the typical 
size of the corona. We assume R = 10R g . Then, on average 
for the 10 observations, the total compactness hot, he ~100, 
and the soft compactness l s ,hc ~ 15 for the hot corona, 
while ltot,wc ^230 and l StWC ~115 for the warm corona. 
For such values of ltot and l s (and thus lh = hot ~ Is the 
heating compactness), no pair equilibrium can be reached 
for the WARM corona giv en its very low temperature (cf. 
iGhisellini fc Haardtlll99l G94 hereafter, and their Fig. 2). 
This is different for the hot corona, which has a high 
enough temperature, given its values of ltot and l s , to be 
dominated by pairs. 

For a pure pair plasma, there is a one-to-one relationship 
between the hard and soft compactnesses and the spectral 
characteristics (energy spectral index a and temperature 
kT) of the Compton spectrum. G94 have mapped the plane 
a-kT into the plane (lh/ls-lh)- We have reported their Fig. 
2 in our Fig. I12[ overplotting the values of lh,hc and l s .h c 
computed for the 10 observations of Mrk 509. 

Astonishingly enough, given the simple assumptions 
of the G94 model (homogeneous spherical corona in a 
homogeneous black body photon field) and our approxi- 
mated estimates of the compactnesses, they agree nicely 
with the theoretical expectations of a pure pair plasma. 
The observed spectral index a ^0.7-0.8 (see Table 1, with 
a = r — 1), the hot corona temperature on the order of 
100 keV, and the variations of these parameters during 
the monitoring (i.e. Aa ~ 0.1 and AfcT ~ 80 keV) are not 
far from the expected ranges. While this comparison is ad- 
mittedly rough (e.g. the observed spectral index a should 
have been steeper between 0.8 and 0.9 to completely agree 
with the expectation), it is interesting to notice that the 
HOT corona behaves spectrally as if it is dominated by pairs. 

In the case of a pair-dominated plasma, the increase in 
the cooling, for a constant heating, is expected to produce 
an increase in the plasma temperature and not a decrease. 
This can be seen in Fig. 1121 Fixing lh but increasing l s 
shifts the pair plasma to a higher temperature state. This 
can be understood as follows (see G94): increasing the cool- 
ing first decreases the plasma temperature. But this implies 
a decrease in the pair creation rate, i.e., a decrease in the 
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Fig. 12. Hard-to-soft compactness ratio Ih/h vs. 
hard compactness l s . T he figure is extracted from 
IGhisellini & Haardtl (|l994h . It shows the one-to-one cor- 
respondence between the spectral characteristics of the 
Compton spectrum of a pure pair plasma (energy spectral 
index a and temperature kT) and the soft and hard 
compactnesses. The black points correspond to the 10 
observations of Mrk 509. 



number of particles. The available energy is thus shared 
among fewer particles and eventually, when pair equilib- 
rium is reached again, the temperature is higher. 

From our fits, the soft compactness l s ,hc of the hot 
corona varies by ^60% during the monitoring, while the 
heating compactness lh,hc varies by ~15%. Thus we are 
roughly in the case of a constant heating with variable 
cooling. Interestingly enough, the hot corona temperature 
deduced from our modeling (marginally) correlates with 
h,hc as expected for a pair-dominated plasma. The linear 
Pearson correlation coefficient is equal to 0.67, which cor- 
responds to a 95% confidence for its significance. 

6.6. Some hint of the behavior on short time scale? 

To investigate the short time scale correlations between 
spectral components in the X-ray band, we estim ated the 
coherence function (e.g. IVaughan k, Nowaklll997t ) between 
the 0.3-1 keV and 1-10 keV energy bands by making use of 
the XMM-Newton data (see right plot of Fig.H"3|). The two 
energy bands display a high level of linear correlation over 
all the sampled frequencies (corresponding to time scales 
between 60 and 0.5 ks). This behavior is at odds with the 
observed lack of correlated variability on longer time scales 
(see Fig. [2] and Sect. [3]). However, this result is consistent 
with the existence of a small and compact hot corona, 
which is able to produce coherent variations on short time 
scales. 

We explored the causal connection between the different 
spectral components by computing time lags as a function 
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Fig. 13. Left: Lag energy spectrum. It is obtained by computing time lags in Fourier frequency space, between the 
energy band 0.3-10 keV and a number of adjacent energy channels, averaging the resulting lags over the frequency interval 
10 -5 -10 -4 Hz. Right: Coherence between 0.3-1 and 1-10 keV energy band s as a function of Fourier frequency once a 
correction for the Poisson noise effects is applied (jVaughan fc NowakL Il997h . 



of energy (see left plot of Fig. [TBI . The time lags are first 
computed in the Fourier frequency domain between the 
entire 0.3-10 keV reference band and adjacent channels, 
with the channel being subtracted from the reference at 
each step, so as to remove the contrib ution from correlated 
Poisson noise (e.g. lZoghbi et al.ll201lD . The derived lags are 
then averaged over the range of frequencies At/ ~ 10~ 5 — 
10 -4 Hz, i.e., the frequency domain not dominated by 
Poisson noise, and plotted as a function of energy. The lag 
spectrum has been shifted so that the zero-lag energy chan- 
nel corresponds to the leading one. Hard positive lags (vari- 
ations in the hard channels being delayed with respect to 
those in the soft channels) are observed over the full energy 
band, the relative amplitude of the measured lags increas- 
ing with the energy separation between the channels. This 
trend is commonly observed in black hole X-ray binaries, 
as well as in some AGNs in the low-frequency range (e.g . 
Nowak et al.l 1 19991: iPapadakis. Nandra. fc Kazanasl [2001; 
Arevalo et all 120061: IZoghbi et all 1201 H lUttlev et al J 12011 
De Marco et al.ll2012D . 

These hard positive lags are generally interpreted by the 
prop agation of mass accretion rate fluctua tions in the disk 
(e.g. lKotov et alj200llArevalo et aT]|2006D and, in order to 
be observed in the soft X-ray bands, the soft X-ray emitting 
region has to be strongly connected somehow to the disk 
accretion flow. This is naturally expected if the soft X-ray 
emitting region is the upper layer of the accretion disk, as 
we propose for the WARM corona. 



6.7. Do we have all the pieces of the puzzle? 

Most of our results are model-dependent. For instance, we 
could add other spectral components (like e.g. a blurred 
reflection or a pure disk black body) to fit the UV and 
soft X-rays. The possibilities are numerous and may give a 
different picture of the inner regions of Mrk 509 than the 
one revealed by our modeling. Our choice was instead to 
choose a small number of spectral components. Our main 
assumption was that the optical-UV emission, observed in 
the optical monitor of XMM-Newton, was the signature of 



an optically thick multi-color accretion disk whose optical- 
UV photons are comptonised in two different media: the 
first one (the warm corona) producing the soft X-ray ex- 
cess and the second one (the hot corona) the high-energy 
continuum emission. In this framework it appears that 

— The warm corona has an amplification ratio equal to 
two i.e. the expected value of an optically thick plasma 
above a passive accretion disk, the disk emission being 
dominated by the reprocessing. In consequence, most of 
the disk emission is due to the heating from the warm 
corona itself. The warm corona could be the warm up- 
per layer of the accretion disk. 

— The hot corona amplification ratio is between five and 
ten and thus suggests a more photon-starved geometry. 
The temperature of the soft seed photons entering and 
cooling the hot corona, on the order of 100 eV, suggests 
a location in the inner region of the accretion flow. 

— A 100 eV temperature is about a factor 5 above the 
expected value for a standard accretion disk around 
a 1.4xlO 8 M0 black hole. Effects of electron scatter- 
ing and of special and general relativity could signifi- 
cantly affect the observed emission and lead to mod- 
ified black body spectra with an "observed" tempera- 
ture higher than the effective temperature. The seed soft 
photons of the hot corona can also be part of the warm 
corona emission. The 100 eV temperature obtained by 
our xspec modeling would be explained by the fact 
that xspec would try to mimic the warm plasma emis- 
sion, which covers an energy range between 3 and ~ 1 
keV, with the multicolor-disk spectral energy distribu- 
tion available in COMPPS for the seed soft-photon field. 

— The HOT corona may be dominated by pairs. 

— The observed flux and spectral variability are mainly 
due to changes in coronae geometries. 



The toy picture presented in Sect. 16.2.31 broadly agrees 
with the non-detection of a relativistically broadened iron 
line in Mrk 509. Indeed, P12 found that the resolved com- 
ponent of the iron line profile (cr=0.22 keV) agrees with an 
origin beyond 40i? g from the black hole. This supports the 
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absence of a neutral optically thick accretion disk down to 
Risco- On the other hand, a weak ionized iron line com- 
ponent is present and can be modeled with a relativistic 
line produced, in an ionized disk, down to ~10 R g from the 
black hole. The presence of this ionized matter agrees with 

a WARM upper layer on the accretion disk; 

It has been recently suggested by iDone etaD (l20Tl 
that the accretion disk in AGNs could evolve radially from 
an optically thick accretion disk beyond a transition radius 
called Rcorona to a two-phase accretion flow at smaller 
radii, formed by an optically thick Comptonized disk and 
an optically thin corona. The latter would play the role 
of our hot corona, while the optically thick Comptonized 
disk will produce the warm Comptonized emission. The 
energetic coupling between the two coronae is taken 
into account in this model. It possesses the different 
components used in our fits (the accretion disk and the 
warm and HOT plasmas). However, it assumes that the 
warm and hot coronae are localized in the same region of 
the accretion flow, below Rcorona-, with the WARM coron a 
embedded in the hot one (see Fig. 5 of lDone et aLl[2uTl . 
This picture does not agree completely with our estimates 
of the amplification ratios; however, it is interesting to 
see that our conclusions are in broad agreement with the 
assumptions proposed by these authors. It will be worth 
testing their model on the Mrk 509 dat£@. This is devoted 
to a future work. 



7. Conclusions 

Mrk 509 was monitored simultaneously with XMM-Newton 
and INTEGRAL for one month in late 2009, with one 
observation every four days for a total of ten observations. 
A series of papers has already been published using these 
data focusing on different aspects of the campaign. Here we 
analyzed the broad-band optical/UV/X-ray /gamma-ray 
continuum of the source during this monitoring. Our main 
effort was to adopt realistic Comptonization models to fit 
the primary continuum, the accretion disk, and the soft 
X-ray excess. 

We obtained a relatively good fit to the data by 
assuming the existence of two clearly different media: a hot 
(kT ~ 100 keV), optically-thin (r ~ 0.5) corona for the 
primary continuum emission, and a warm (kT ~ 1 keV), 
optically-thick (r ~ 10-20) plasma for the optical-UV to 
soft X-ray emission. 

These two media have very different geometries: the 
warm corona covers a large part of the accretion disk, the 
disk radiation being dominated by the reprocessing of the 
warm corona emission. The warm corona could be the 
warm upper layer of the accretion disk. The temperature 
of the soft-photon field cooling it is on the order of 3 eV. 

In contrast, the HOT corona has a more photon-starved 
geometry and can be the inner part of the accretion flow. 
The transition radius Ri„ between the outer disk/ warm 
plasma and the HOT inner corona is estimated to be around 
10-20 R g . The variation in M but also in Ri n would be 
at the origin of the variation in the coronae luminosities, 

4 Their model is called OPTXAGN, and is publicly available as 
a local model for XSPEC 



the variation in Ri n also producing the variation in the 
hot corona amplification ratio. The spectral and flux 
variation in the HOT corona are also consistent with a 
pair-dominated plasma. 

These different conclusions may depend at some level on 
the choice of our modelling. For example, we do not include 
blurred ionized reflection, while part of the soft X-ray excess 
could be produced by this component. The presence of an 
(admittedly weak) ionized iron line, which agrees with a 
relativistic profile, by P12 would also support the presence 
of such reflection. On the typical time scale (~ day) of 
this campaign, however, it does not seem to play a major 
role and Comptonization in a warm plasma, covering the 
accretion disk, appears the most natural explanation for 
the soft X-ray excess in Mrk 509. 
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Appendix A: Optical depth estimates 

The NTHCOMP model used to fit the optical-UV-soft X-ray 
emission does not directly provide the value of the optical 
depth of the plasma for which it models the emission. 
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However, the plasma optical depth has more physical 
meaning than the photon index of the emitted spectrum, 
so it is an important parameter that we need to estimate. 
Generally, one compares the nthcomp spectral shape, for 
a given set of photon index and temperature, to another 
model that provides the optical depth. While commonly 
used, the procedure is generally not precisely discussed. 
The accuracy of this comparison is important, because 
inaccurate values for the temperature and optical depth 
will give a different interpretation for the heating/cooling 
ratios (see Sect. I6.1[) . 



We compared the results obtained with nthcomp 
with those from three other xspec Comptonization 
models: COMPST, COMPT T, and eqpair. COMPST 
(jSunvaev k. Titarchuld . Il980t) only has the corona temper- 
ature and optical depth as free parameters (no parameter 
for the soft photon temperature), but the range of validity 
is relatively wide, including low temperature and large 
optic al depth as expected in our case. COMPTT (|Titarc huk. 
1994) is a commonly used Comptonization code that 
provides the optical depth and the temperature. Like 
COMPST it assumes a black body distribution for the soft 
photon field. This is the model used by Mil. eqpair is 
a more complete Comptonization model that carefully 
takes most of the microphysics needed to compute re alistic 
Comptonization spectra into account (jCoppil . 119911) . The 
soft photon distribution can be cither a simple black 
body or a multicolor disk spectrum. It can treat hybrid 
thermal/non-thermal particle distributions. Its natural 
parameters are the optical depth and the heating/cooling 
ratio. The temperature is an output of the code. To obtain 
the temperature in xspec, the chatter command has to 
be larger than 15. However, for most of these models, the 
values expected for the temperature and optical depth of 
the WARM corona are outside the nominal range of validity 
indicated in xspec. The temperature has to be between 
2 and 500 keV for COMPTT, between 1 keV and 1 MeV 
for nthcomp, and the optical depth should be below 10 
for eqpair. It is thus interesting to see how these models 
compare. 



We report in Figs. lA.lb . b, and c the contours of 
the temperature and optical depth for the three models 
COMPST, COMPTT, and eqpair in the T-kT nt h parameter 
space of NTHCOMP. We restrict our analysis to low temper- 
atures and steep spectra as expected for the WARM corona. 
While the contours of COMPST and EQPAIR look relatively 
similar, the optical depth estimated with COMPTT is 
a factor two smaller. The corresponding optical depth 
estimates for the ten different observations of Mrk 509 
discussed in this paper are plotted in Fig. lA.lb !. Possible 
explanations of this difference can be found in the different 
model assumptions, especially concerning the geometry of 
the disk-corona systems. It is also worth noting that the 
best fit temperature values of the warm corona obtained 
by Mil are significantly lower (by a factor 4) than ours. 
In this case, however, the main reason for this difference 
is certainly in the use of a different model for the primary 
continuum (Mil use a power law). Interestingly enough, 
while the exact values of the optical depth clearly depend 
on the model used (cf. Fig. lA.lt i). its variability through- 
out the monitoring looks very similar and is apparently 



independent of the model used. 



Appendix B: Estimation of the L s and L tot 

To estimate the seed soft photon luminosity entering and 
cooling a corona of optical depth r, we use the conservation 
of the number of photons, which characterizes the Compton 
process. From the observed comptonized spectrum, we de- 
duce an observed photon rate n b s . By conservation, this 
photon rate is equal to the sum of the seed photons cross- 
ing the corona without being comptonized rt s o and those 
comptonized in the corona and emitted upward (in direc- 
tion to the observer) n CjUp . At first order, 



n s ,o 



n s e 
n s (l 



(B.l) 
(B.2) 



the factor 1/2 of the last equation assumes the Compton 
scattering process to be isotropic with half of the comp- 
tonized photons being emitted upward, and the over half 
being emitted downward. The Compton process is certainly 

he case of an anisotro pic seed 
rdt fc Maraschil [19911 lHaardtl 
but the 

effect on the photon rate is relatively small. 



not isotropic, especially in the case of an an isotro pic seed 

-9911: 

119931 : iStern et al.1 [l995t iHenri fc Petruccil 11991 7 



soft photon field (e.g. iHaarc 



Then we have 



n s . 



which gives 



l 



(B.3) 



(B.4) 



For an optically thin corona (r < 1), n, ~ n b s (very few 
soft photons are comptonized). This is rougly the case of 
the hot corona. On the other hand, for r > 1, n, ~ 2n t, s , 
which corresponds to the case of the WARM corona. Once 
the photon rate n s of the multicolor accretion disk is 
known, as is its temperature (deduced from the fit), we 
have access to the corresponding luminosity L s . 

The corona total power L to t is the sum of L s and the 
heating power L\ x liberated in the corona to comptonize the 
soft photon field: 

L tot = L h + L s . (B.5) 

Here, Lh and L s can be divided into upper and lower parts: 

Lh = Lj-^u + Lh,d = ZLh^u (B-6) 
L x = L s 



, L s (l-e- T ) L s (l-e- T ) 



(B.7) 
(B.8) 



The first equation assumes an isotropic Compton process, 
and the second equation is obtained with the same reason- 
ing as in Eq. |BT] and IIOI 



Then the observed luminosity is simply 



(B.9) 
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Fig. A.l. Contours of the temperature (dashed lines) and optical depth (solid lines) in the T-kT nt h parameter space of 
NTHCOMP for COMPST (top left), COMPTT (top right) and eqpair (bottom left) respectively. Bottom right: Light 
curves of the optical depth estimates during the whole campaign for the three models: diamonds: COMPST, squares: 
eqpair; triangles:COMPTT 



where we neglect any extra emission that does not come 
from the corona. From Eqs. IB. 61 IB. 81 and IB.91 we can de- 
duce Lh 



L h = 2L obs - L s (l + e r ), 
and finally, from Eq. IB. 51 
Ltot = 2L h s — L s e T . 



(B.10) 



(B.11) 



For an optically thin corona (t <C 1), L to t — 2Lobs — L s . 
This is roughly the case for the hot corona. On the other 
hand, for r > 1, L to t — 2L (, S , which corresponds to the 
case of the warm corona. 

The radiative equilibrium of the soft phase (which pro- 
duces the soft photon field) also implies 



L s ,d + Lh t d + L s ,intr 

2 2 s,intr 



(B.12) 
(B.13) 



where L s j ntr is the intrinsic emission of the soft phase. In 
consequence, the heating/cooling ratio is equal to 



Lf. L s 



(B.14) 



In the case of an optically thin (r < 1) corona above a 
passive disk (L s .j ntr = 0), we found the well kno wn result 
Lh/L s — 2 (see e.g. Haardt fc Maraschil 1 199 ll and their 
Eq. 3b with f=0 and no disk albedo) and consequently 
L to t/L s = 3. On the other hand, for an optically thick 
corona (still above a passive disk), we find Lh/L s = 1 and 
Ltot/ L s = 2 
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